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We report improved measurements of B decays with an n meson in the final state using 492 fb _1 of 
data collected by the Belle detector at the KEKB e + e~ collider. We observe the decays B — > lyx 
and B — > tjK^ and measure the branching fractions B(B ± — > T)ir ) = (4.2±0.4(stat) ±0.2(sys)) x 
1(T 6 and B(B ± -» TfK*) = (1.9 ± 0.3(stat)+°;*(sys)) x 10" 6 . The corresponding CP-violating 
asymmetries are measured to be —0.23 ± 0.09(stat) ± 0.02(sys) for jpr^ and —0.39 ± 0.16(stat) ± 
0.03(sys) for rjK ± . We also search for B° — > r]K° decays and set an upper limit of 1.9 x 10~ 6 at the 
90% confidence level. 
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Charmless B decays provide a rich sample to under- 
stand B decay dynamics and to search for CP viola- 
tion. The decay B — > r]K proceeds through a b — *■ s 
penguin process and a b — ► u tree transition. Inter- 
ference from the two penguin processes, b — > sss and 
b — ► uus, and the known rj — rj' mixing are expected to 
enhance the B — ► vj K branching fraction but suppress 
The situation is reversed for the r]K* and 
rj'K* modes since the r] and K* mesons are in a rela- 
tive p-wave rather than an s-wave state. Experimental 
results @, d, 0, [B| have confirmed this picture but more 
precise measurements of r/K and rj'K* are needed for a 
quantitative understanding of the underlying dynamics. 
Moreover, the penguin amplitude of r]K may interfere 
with the b — > u amplitude, resulting in a large direct CP 
asymmetry (Acp) @. Theoretical expectations for con- 
tributions from other mechanisms @, d, Q also suggest a 
large Acp although the sign could be positive or nega- 
tive. Our earlier measurements with limited statistics [3( 
indicated a large negative Ac p central value for r]K ± . 

The study of B° — > r\K is of particular interest be- 
cause this decay is a CP eigenstate and could be used 
for time dependent CP measurements. The dominant 
process in B — * 7777 decays is the (external) b — > u tree 
while a suppressed b — > d penguin process may also 
contribute. It has been argued [7|, llfj that the direct 
CP violating asymmetry could be large in the 777r ± and 
rfn^ modes, whose branching fractions are expected to 
be around (2 - 5) x 1(T 6 @,|lo|. 

In this paper, we report improved measurements of 
branching fractions and partial rate asymmetries for 
B — > rjh decays, where ft, is a charged or neutral K meson 



or a charged it meson. The partial rate asymmetry for 
charged B decays is defined to be: 

_ N(B- ^ V h-) - N(B+ ^ V h+) 
CP N(B~ — > r]h~) + N(B+ — > r]h + ) ' U 

where N(B~ — > r]h~) is the yield obtained for the 
B~ — > r)h~ decay and N(B + — > 77 h + ) denotes that of 
the charge-conjugate mode. The data sample consists 
of 535 million BB pairs (492 ftr 1 ) collected with the 
Belle detector at the KEKB e + e~ asymmetric-energy 
(3.5 on 8 GeV) collider [ll| operating at the T(4S) res- 
onance. Throughout this paper, the inclusion of the 
charge-conjugate decay mode is implied unless otherwise 
stated. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect K\ mesons and to identify muons (KLM) . The 
detector is described in detail elsewhere [l2|. In August 
2003, the three-layer SVD was replaced by a four-layer 
device with greater radiation tolerance 1 1 31 ] - The data 
sample used in this analysis consists of 140 fb _1 of data 
with the old SVD (Set I) and 352 fb _1 with the new one 
(Set II). 

The event selection and B candidate reconstruction are 
similar to those documented in our previous publication 
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Two 77 decay channels are considered in this analy- 
sis: 77 — ♦ 77 (t7 77 ) and 77 — ► 7r + 7r~7r (773-^). We require 
photons from the 77 and 7r° candidates to have laboratory 
energies (Ej) above 50 MeV. In the ?7 77 reconstruction, 

the photon energy asymmetry, ^rrg-— , is required to 
be less than 0.9 to reduce the large combinatorial back- 
ground from soft photons. Neither photon from ry 77 is 
allowed to pair with any other photon having £" 7 > 100 
MeV to form a ir° candidate. Candidate tt° mesons are 
selected by requiring the two-photon invariant mass to be 
in a mass window between 115 MeV/c 2 and 152 MeV/c 2 . 
The momentum vector of each photon is then readjusted 
to constrain the mass of the photon pair to the nominal 
7T° mass. 

Candidate 773^ mesons are reconstructed by combining 
7T° candidates with at least 250 MeV/c laboratory mo- 
mentum with a pair of oppositely charged tracks that 
originate from the interaction point (IP). We impose the 
following requirements on the invariant mass of the 77 
candidates in both data sets: 516 MeV/c 2 < M 77 < 569 
MeV/c 2 for ?7 77 and 539 MeV/c 2 < M 37r < 556 MeV/c 2 
for 773^. After the selection of each candidate, the 77 mass 
constraint is implemented by readjusting the momentum 
vectors of the daughter particles. 

Charged tracks arc required to come from the IP. 
Charged kaons and pions, which are combined with 77 
mesons to form B candidates, are identified using a K{ir) 
likelihood LK(L n ) obtained by combining information 
from the CDC (dE/dx), the TOF and the ACC. Dis- 
crimination between kaons and pions is achieved through 
a requirement on the likelihood ratio Lk/{L-„ + Lk)- 
Charged tracks with likelihood ratios greater than 0.6 are 
regarded as kaons, and less than 0.4 as pions. Charged 
tracks that are positively identified as electrons or muons 
are rejected. The K/tt identification efficiencies (PID) 
and misidentification rates are determined from a sam- 
ple of D* + — > D°TT + ,D a — > K~ir + decays with kaons 
and pions in the same kinematic region of two-body B 
decays. The kaon (pion) identification efficiency is 83% 
(90%) and 6.4% (11.7%) of pions (kaons) will be misiden- 
tified as kaons (pions). The systematic error of the K/tt 
selection is about 1.3% for pions and 1.5% for kaons, re- 
spectively. 

Kg candidates are reconstructed from pairs of 
oppositely-charged tracks with an invariant mass (M n7r ) 
between 480 MeV/c 2 and 516 MeV/c 2 . Each candidate 
must have a displaced vertex with a flight direction con- 
sistent with that of a ifly-meson originating from the IP. 

Candidate B mesons are identified using the beam- 
energy constrained mass, Mb c = \J E'£ cllm — P B , and the 
energy difference, AE = E B - E heam , where E hca-m is 
the run-dependent beam energy in the T(45) rest frame 
determined from B — > D^n events, and Pb and Eb are 
the momentum and energy, respectively of the B candi- 
date in the T(45) rest frame. The resolutions in Mb c and 



AE are about 3 MeV/c 2 and 20-30 MeV, respectively. 
Events with M hc > 5.2 GeV/c 2 and \AE\ < 0.3 GcV are 
selected for the analysis. 

The dominant background comes from the e + e~ — > qq 
continuum, where q = u, d, s or c. To distinguish sig- 
nal from the jet-like continuum background, event shape 
variables and B flavor tagging information are employed. 
We combine the correlated shape variables into a Fisher 
discriminant 14| and then compute a likelihood that is 
the product of probabilities based on this discriminant 
and cosOb, where 9b is the angle between the B flight 
direction and the beam direction in the T(45 l ) rest frame. 
A likelihood ratio, 1Z = C S /(C S + C g g), is formed from 
signal (C s ) and background (C q q) likelihoods, obtained 
from Monte Carlo simulation (MC) and from data with 
Mbc < 5.26 GeV/c 2 , respectively. Signal MC events for 
the charged B modes are generated with the PHOTOS 
[TBI simulation package to take into account final state 
radiation. Additional background discrimination is pro- 
vided by B flavor tagging. Events that contain a lepton 
(such as those used in high quality tagging) are more 
likely to be BB events so a looser 1Z requirement is ap- 
plied. The standard Belle B tagging package [lj| pro- 
vides two outputs: a discrete variable (q) indicating the 
tagged side flavor and a dilution factor (r) ranging from 
zero for no flavor information to unity for unambiguous 
flavor assignment. Since the charged B modes are fla- 
vor specific, the wrong flavor tagged events are likely to 
be background and a tight 1Z requirement can be ap- 
plied. We divide the data into six sub-samples based on 
the q and r information for the charged modes and the 
r value only for the neutral mode. Continuum suppres- 
sion is achieved by applying a mode-dependent require- 
ment on 1Z for events in each sub-sample that maximizes 

N™ p /yjNs Kp + N™ p , where A s ox p is the number of signal 

events expected from MC and iV°| p denotes the number 
of background events estimated from data. After apply- 
ing the 1Z requirements, we select one candidate per event 
based on the best 1Z. The fraction of events with multi- 
ple candidates are ~ 1% for the 77 mode and ~2-3% for 

the 1T + TT~ 



mode. 



Using a large MC sample, all other backgrounds are 
found to be negligible except for r]K + {rjTT + ) reflecting 
into the r/ir + (r]K + ) sample, due to K + — ► h + (tt + — ► K + ) 
misidentification, and the feed-down from charmless B 
decays, predominantly B — ► 77if*(892) and B — ► ?]p(770). 
We include the reflection and charmless components in 
the fit used to extract the signal. 

The signal yields and partial rate asymmetries are ob- 
tained using an extended unbinncd maximum-likelihood 
(ML) fit with input variables Mb c and AE. The likcli- 
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hood is defined as: 

C = e-^i N * x nCC^j) and ( 2 ) 

i 3 

V) = ^[l-q'-AcP^PjiMLAE*), (3) 

where i is the identifier of the z-th event and Nj is the 
number of events for the category j, which corresponds 
to either signal, qq continuum, the reflection due to K-tt 
misidentification, or background from other charmless B 
decays. Pj{M\, c , AE) is the two-dimensional probability 
density function (PDF) in Mb c and AE, and q indicates 
the B meson flavor, B + (q = +1) or B~(q = —1). For the 
neutral B mode, Vj in Eq. 2 is simply Pj{M^ c , AE 1 ) and 
there is no component from charged particle misidentifi- 
cation. 

In Ref. we reported that in both data sets the PID 
efficiency is slightly different for positively and negatively 
charged particles. Therefore, the raw asymmetry defined 
in Eq. 1 must be corrected. This efficiency difference 
results in an A C p bias of -0.005 (+0.005) for rjir (r)K). 
The bias is subtracted from the raw asymmetry. 

The PDFs for the signal, the reflection background 
and the charmless feed-down are modeled with two- 
dimensional Mbc-AE smooth functions obtained using 
large MC samples. The signal peak positions and resolu- 
tions in Mbc and AE are adjusted according to the data- 
MC differences using large control samples of B — > Dir 
and D° — > K + it~tt° decays. The continuum background 
in AE is described by a first- or second-order polynomial, 
while the Mb c distribution is parameterized by an AR- 
GUS function, f(x) = x\[\ — x 2 exp [— £(1 — x 2 )], where 
x is Mbc/-Ebcam \2M ■ The continuum PDF is thus formed 
by the product of an ARGUS function and a polynomial, 
where £ and the coefficients of the polynomial arc free 
parameters. 

The partial rate asymmetries of the charmless B back- 
grounds are fixed to zero in the fit while the Acp and 
normalizations of the reflection components are fixed to 
expectations based on the B + — > r/K + and B + — > rjir + 
partial rate asymmetries and branching fractions, as well 
as K + «-» 7T + fake rates. The reflection yield and Aqp 
are first input with the assumed values and are then re- 
calculated according to our measured results. 

Table fl] shows the measured branching fractions for 
each decay mode as well as other quantities associated 
with the measurements. The efficiency for each mode is 
determined using MC simulation and corrected for the 
data-MC discrepancy obtained from the control sample 
studies. In addition to the particle identification perfor- 
mance discrepancy, our MC slightly overestimates the ef- 
ficiency for detecting low momentum 7r°s, which results 
in a 3.1% correction for the mode. The combined 
branching fraction for the two datasets is computed as 
the sum of the yield divided by its efficiency in each set 
divided by the number of B mesons, while the partial 



rate asymmetry for the charged mode is computed using 
the sum of the yield divided by its efficiency in each set in 
Eq. 1. The combined branching fraction and partial rate 
asymmetry of the two rj decay modes are obtained from 
the weighted average assuming the errors are Gaussian. 
The number of B + B~ and B°B° pairs are assumed to 
be equal. Figure Q] shows the M^ c and AE projections 
after requiring events to satisfy —0.10 GeV < AE < 0.08 
GeV and M bc > 5.27 GeV/c 2 , respectively. 
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1: Mbc and AE projections for (a,b) B — ► rpr , (c,d) 
B ± — > rjK ± , and (e,f) B° — > r/K decays with the r]~ ll and 
773^ modes combined. Open histograms are data, solid curves 
are the fit functions, dashed lines show the continuum contri- 
butions and shaded histograms are the feed-down component 
from charmless B decays. The small contributions around 
Afbc = 5.28 GeV/c 2 and AE = ±0.05 GeV in (a)-(d) are the 
reflection backgrounds from B — > T]K ± and B — > 77^*. 

Systematic uncertainties due to the signal PDFs used 
in the fit are estimated by performing the fit after vary- 
ing the signal peak positions and resolutions by one stan- 
dard deviation (a). We also examine the changes in yield 
and Acp when the requirement of no asymmetry for the 
charmless background is removed. In B ± — > ?77r ± , the 
reflection yields are estimated to be 9.4 ± 3.1 events for 
the ?7 77 mode and 3.6 ± 1.9 for 773^ while in B^ — * rjK^ , 
the reflection yields are 13.9±3.7 for ?7 77 and 4.6±2.1 for 
773^ . The reflection yields and their Acp values are varied 
by one standard deviation in the fit to obtain the corre- 
sponding systematic uncertainties. The quadratic sum of 
the deviations from the central value gives the systematic 
uncertainty in the fit. A statistical significance is calcu- 
lated as S = \J — 21n£o (— 21n£ max ), where — 21n£o 
is for zero signal yield and — 21n£ max is for the best-fit 
value. The final significance including systematic uncer- 
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TABLE I: Detection efficiency (e) including sub-decay branching fraction [IjJ, yield, significance (Sig.), measured branching 
fraction (£>), the 90% C.L. upper limit (UL) and Acp for the B — > r\h decays. The first errors in columns 3, 5 and 7 are 
statistical and the second errors are systematic. 



Mode e(%) Yield Sig. g(l(T H ) UL(KT K ) A C p 



B ± -> 7 7 7T ± 






14.7 


4.2 ±0.4 ±0.2 




-0.23 ±0.09 ±0.02 


77 77 7T ± 


8.3 


183 ± 20 


11.9 


4.1±g;| ±0.2 




-0.11 ±0.11 ±0.01 




3.1 


73111 


8.7 


4.4l^ ±0.3 




-0.52 ±0.16 ±0.02 


B ± -> r/A* 






8.1 


1.9±0.3loi 




-0.39 ±0.16 ±0.03 


r/ 1 ~ / K ± 


7.3 




6.4 


1.9±°;*±0.1 




-0.30 ±0.19 ±0.02 


r?3^A± 


2.7 


29 ±8 


4.8 


2.0±g;|±0.2 




n c;c+0.27+0.05 
U.OO_ .28-0.04 


S° -» 77A 






2.9 


1.1 ±0.4 ±0.1 


< 1.9 




77 77 A° 


2.6 


16 ±8 


2.6 


l.l±g;S ± 0.1 


< 2.2 






1.0 




1.2 


0.9±g;| ±0.1 


< 2.4 





tainty is taken as S = S — \/Y>(S — Si) 2 , where S is 
the statistical significance for the fit and Si is the signif- 
icance obtained for each systematic check with the value 
smaller than S . 

The possible detector bias due to the tracking accep- 
tance for Acp(B ± — > Tj7T^~ ) and Acp{B ± — > ^A^) is 
evaluated using the Acp value of the continuum com- 
ponent. No obvious bias is observed and we use the 
statistical error of the ?7 77 and 773^ modes combined as 
the systematic error. The bias error of 0.01 is added in 
quadrature with the fit systematic error to give the final 
systematic uncertainty in Table I. Figures [2] and [3] show 
the Mbc and AE projections for the B + and B~ samples. 
In both the rjir^ and r\K^- modes, we observe larger B + 
signals. 

The systematic error of the efficiency arises from the 
1Z requirement, tracking efficiency, particle identification, 
Kg reconstruction, tt and ?y 77 reconstruction, and r? 77 
and 773^ branching fractions. The performance of the 
1Z requirement is studied by checking the data-MC ef- 
ficiency ratio using the B + — > D°tt + control sample. 
The systematic errors on the charged track reconstruc- 
tion are estimated to be 1% per track using partially 
reconstructed D* events. The ir° and ?7 77 reconstruction 
efficiency is verified by comparing the tt decay angu- 
lar distribution with the MC prediction, and by mea- 
suring the ratio of the branching fractions for the two 
D decay channels D° — ■> 7\ + 7r~ and D° — > A' + 7r~7r°. 
The Kg reconstruction is verified by comparing the ra- 



tio of D+ 



K-ir+ir+ yields. The 



uncertainties in the ?/ 77 and 773^ branching fractions are 
taken from Ref. Table [TT] summarizes the systematic 
uncertainties, including the error on the number of BB 
events. The systematic error that arises from how well 
PHOTOS describes final state radiation is found to be 
negligible [2(j. The final systematic error for the com- 
bined branching fraction is obtained by assuming that 
the systematic errors for the sub-decay modes are 100% 
correlated. 

We observe an excess of B° — > r/K events but the 
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FIG. 2: Mbc and AE projections for (left) B~ — ► 7771-" and 
(right) B + — » rj-K + with the ?7 77 and 773^ modes combined. 
Open histograms are data, solid curves are the fit functions, 
dashed lines show the continuum contributions and shaded 
histograms are the contributions from charmless B decays. 
The small contributions near M bc = 5.28 GeV/c 2 and AE = 
—0.05 GeV are the backgrounds from misidentified B — * 
r]K ± (reflections). 



significance is slightly less than 3. Therefore, an upper 
limit on the branching fraction at 90% confidence level 
is provided. To calculate this limit, we find the yield for 
which 90% of the area of the likelihood function lies at 
lower values. We divide the yield by the reconstruction 
efficiency reduced by la of its uncertainty, which is the 
quadratic sum of the errors given in Rows 2-8 of Table 
2. The result is then inflated by the lcr uncertainty due 
to the parameters fixed in the fit (1st row of Table 2) 
to obtain the upper limit including all systematic uncer- 
tainties. 
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TABLE II: The systematic uncertainties for the B — » rjh branching fractions, given in %. The fit systematic errors include the 
uncertainties due to the signal PDFs, the yields of the reflection backgrounds and the partial rate asymmetries of the charmless 
B and reflection backgrounds. 



Sources 


rj- n ir ± 


»737r 7T* 


■q 11 K ± 




r ?77 iC u 


773* K u 


Fit 


+2.6 

—3.0 


+ 3.0 
— 3.6 


+6.1 

—4.3 


+6.2 
—6.4 


+6.0 
—6.2 


±6.7 


1Z requirement 


±1.2 


±1.2 


±1.2 


±1.6 


±1.4 


±1.4 


Tracking 


±1.0 


±3.0 


±1.0 


±3.0 




±2.0 


PID 


±1.3 


±1.3 


±1.5 


±1.5 






Kg reconstruction 










±4.9 


±4.9 


77 reconstruction 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


B(r) -> 77) 


±0.7 




±0.7 




±0.7 




B(j] -> TT+TT^TT ) 




±1.8 




±1.8 




±1.8 


N B 


±1.3 


±1.3 


±1.3 


±1.3 


±1.3 


±1.3 


Sum 


+5.4 
-5.6 


+6.b 

-6.6 


+7.7 
-6.4 


+8.5 

-8.7 


+8.y 

-9.1 


±9.8 




> 

CD 



O 
CN 
l_ 

CD 
Q_ 

CO 
CD 



LJ 



5.3 5.25 

M bc (GeV/c 2 ) 



5.3 



: 




1 1 - 

7]\\~ 1 


1 


7]K + 1 









-0.25 



0.25-0.25 

AE (GeV) 



0.25 



FIG. 3: M bc and AE projections for (left) B~ -> r]K' and 
(right) B + -> rjK + with the ?7 77 and 773^ modes combined. 
All symbols are the same as in Fig. 2. The small contribu- 
tions near M bc = 5.28 GeV/c 2 and AE = 0.05 GeV are the 
backgrounds from misidentified B — > ^tt* (reflections). 



the BaBar collaboration are also negative, slightly more 
than lcr away from zero for each mode. Larger data 
samples are needed to verify these large CP asymme- 
tries. Finally, we find a hint of an r]K° signal with 
B{B° -> r/K ) = (1.1 ±0.4 ±0.1) x 10" 6 . Since the mea- 
surement is not significant, we provide an upper limit at 
the 90% confidence level of 1.9 x 10~ 6 . A similar hint was 
also observed by the BaBar collaboration with a central 
value of (1.5 ± 0.7 ± 0.1) x 10~ 6 . The combined aver- 
age, (1.2 ± 0.4) x 10~ 6 , shows 3.4er evidence for the CP 
eigenstate decay B° — > 77 K°. 
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In summary, we have observed B — > r/ir and B — > 
r/K^ decays; their branching fractions are measured to 
be (4.2 ± 0.4 ± 0.2) x 10" 6 and (1.9 ± 0.3^) x 10~ 6 , 
respectively. These results are consistent with our previ- 
ously published measurements Q with statistical errors 
reduced by more than 40%. Compared with the ear- 
lier BaBar results our measurements are more pre- 
cise despite a 1.8cr lower branching fraction on B ± — > 
•qK^ . The CP-violating asymmetries are measured to be 
A CP {B ± -> 7 7 7r ± ) = -0.23±0.09±0.02 and A C p{B ± -> 
rjK^) = -0.39 ± 0.16 ± 0.03, which are 2.5cr and 2.4cr 
away from zero, respectively. It is interesting to note 
that the Acp values for these two modes obtained by 
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